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We consider the feasibility of performing quantum lògic operations based on stimulated Raman 
transitions in trapped Calcium ions. This technique avoids many of the technical difiiculties involved 
with làser stabilisation, and only three làser wavelengths are required, none of which need have 
, particularly stringent requirements on their bandwidths. The possible problems with experimental 

realisations are discussed in detail. 
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INTRODUCTION 



Quantum information processing is a new paradigm in computation, and the possibility of realizing it via trapped 
ions is currently a subject of active research in a number of groups worldwide. The bàsic idea of trapped ion quantum 
' computation is to use two infernal quantum states of a trapped ion to form a "qubit" for the storage of information, and 
t~ — , to employ the strongly coupled quantum oscillations of the ions to perform the conditional dynamics ( or "quantum 
lògic gates" ) required to execute quantum algorithms [ ï], ^) . The quantum dynamics of both the infernal and external 
degrees of freedom of the ions are controlled by pulsed làser beams. Recent reviews of the experimental issues involved 
are given in refs. |3[ |J. 

There are two different ways of performing the laser-induced transitions required for quantum computation. Con- 
ceptually the most simple is to use the Rabi oscillations between two infernal states of the ion. If the ion interacts 
with a monochromatic field (i.e. a single narrow band làser beam) precisely tuned to the transition frcquency for the 
two levels, the population will oscillate back and forth between them ||. A long-lived two level system is required to 
make a qubit suitable for quantum computation, and so atòmic levels with dipole-forbidden transitions are the most 
suitable. An example of such a dipole forbidden transition is the 729 nm 4 2 Si/2 to 3 2 D 5 / 2 transition in Ca + (fig.1.1), 
which has a natural lifetime of about 1.1 seconds ||. We will refer to this method as the "single làser" technique. 
^ \ The theoretical background is described in some detail in ref. || . 

Alternatively, qubits can be manipulated using stimulated Raman transitions. The advantages of using stimulated 
Raman transitions, instead of single-photon transitions, are now well known. Two làsers, traditionally named the 
• i-H . "pump" and the "Stokes" beams, are tuned so that population from one level is excited to some intermediate virtual 
■ level by one làser, and then immediately brought back down from the virtual level to a different atòmic state by 
the second làser. Following the first demonstration of these transitions between atòmic ground state s H , Raman 
transitions were proposed for trapped ion experiments Jïp[ , then used to cool trapped neutral àtoms flïlf . Raman 
transitions have been used in ion traps in a series of experiments by the NIST group |Ï2| ] . In this paper we discuss 
in detail the use of such stimulated Raman transitions to perform quantum lògic operations in ion-trap quantum 
computers based on Ca + , or Ca + -like ions. In particular, we focus on two issues: options for realising quantum lògic 
operations; and measurement of the final state, i.e. the "read out". The latter issue has been investigated in 40 Ca + 
|Q , however that method does not realise an an ideal von Neumann measurement (projection onto an orthonormal 
basis) . 

The principal advantages of the Raman scheme in comparison with the single làser scheme are the much reduced 
optical stability requirements and a reduction in the rate of spontaneous emission JtJ. Only high relative frequency 
stability of the pump and Stokes beams is required, which is easily achieved with current technologies. The two beams 
can be generated by dividing a parent beam into two with a beam-splitter and modulating one of the resultant beams 
using an electro-optic modulator (EOM) driven by a stable radio-frequency (RF) source. Regarding spontaneous 
emission, both the lower and the upper qubit states |0) and |1) are sublevels of the ground state, and so the radiative 
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lifetime of |1) is extremely large (~ 10 6 years or more). As all logical operations must be carried out before spontaneous 
emission of the |1) state can occur, this significantly relaxes the operational constraints. The error rate per CNOT 
operation has been calculated to be two to three orders of magnitude smaller in Raman schemes than in single làser 
transition schemes ||. 

The choice of ion species is driven principally by the availability of reliable làsers at the requisite wavelengths. 
Figure 1.1 shows the level scheme for 40 Ca + , when the 4 2 Si/2 ground state level is Zeeman-split into sublevels by 
a modest magnètic field (the splitting is 2tt x 2.8MHz per Gauss ||). The ground state sublevels of each ion are 
assigned to be the two levels of each qubit (i.e. the logical |0) and |1) states). The two Raman làsers would then 
be used to perform transitions between these two states, via a virtual level. For singly ionised calcium (Ca + ) the 
bàsic wavelength ranges required are 393-397 nm, 729-732 nm and 849-866 nm, which can be realised relativcly casily 
with existing commercially available systems, and (for the 393-397 nm range) a doubling system. Although Ca + is 
attractive in terms of wavelength, the most abundant isotopes of Ca + have zero nuclear spin and thus no hypcrfinc 
splitting. This is unfortunate because the hyperfme structure of the ground state of other species (e.g. 9 Be + , 137 Ba + , 
25 Mg + ) provides both a large number of additional quantum states and the frequency splitting required to perform 
quantum lògic cfhciently using Raman transitions. However, as discussed in section 3, there are a number of options 
available for by-passing this problem. 

The paper is organised as follows. In section 2 we present a concise review of the theoretical description of stimulated 
Raman transitions as applied to trapped àtoms or ions, for which both infernal and external degrees of freedom must 
be considered. Section 3 explains how Raman transitions can be used to perform all of the tasks required for quantum 
computation. Section 4 discusses the experimental issues we anticipate in carrying out this program in the laboratory, 
and section 5 states our conclusions. 



II. THEORY OF STIMULATED RAMAN TRANSITIONS 
A. Basic notions and notation 

Theoretical treatment of stimulated Raman processes in àtoms has been considered by many authors (see, for 
example, the comprehensive works by Schenzle and Brewer flÏ4[ and by Shore For the convenience of the reader, 
we recapitulate the relevant theory here. 

The quantum states of a trapped ion consist of a tensor product of the wavefunction describing the infernal states 
of the ion and of the wavefunction describing the external vibrational state of the ion in the trap. The infernal states 
are specified by a set of quantum numbers; we will denote these states by the short-hand notation I7, J, toj), where 
J, mj are the usual angular momentum quantum numbers and 7 lumps together all of the other quantum numbers 
which specify the state. As discussed below, Calcium has only one very rare stable isotope with non-zero nuclear 
spin. Thus we will assume that we are dealing with spin zero nuclei, so there is no hyperfme structure to consider, 
and Jmj are "good" quantum numbers. We shall divide these infernal states of the ion into "upper" and "lower" 
manifolds. What we have in mind for these manifolds are the various magnètic sublevels of some atòmic level: For 
example, the lower manifold could be the two sub-levels of the 4 2 S;iy 2 ground state, whilc the upper manifold could be 
the two sublevels of the 4 2 P 1/ / 2 or the 4 2 P 3 / 2 excited levels, although we shall develop the theory for arbitrary upper 
and lower manifolds. Further, we will assume that transitions between any state in the upper manifold and any state 
in the lower manifold are dipole allowed, and that transitions between any pair of states in the upper manifold or 
any pair of states in the lower manifold are forbidden. Generalizing the theory to consider the more exòtic case when 
transitions between the upper and lower manifolds are dipole forbidden (e.g. when the upper manifold is the is 3 2 D 3 / 2 
level) is reasonably easy. This situation may have eventual application to ion trap quantum computers, because the 
performance of such devices will be limited by spontaneous emission from the upper manifold, as has been pointed 
out by Plenio and Knight |^5| ; thus if a long-lived upper manifold can be used, in principle higher performance could 
be obtained. 

We will assume that the ions are confined in a highly anisotropic trap, with very strong confining potentials in 
two directions and a comparatively weak potential in the third, axial direction, which we will refer to as the x-axis. 
Provided that the the trapping potentials have a sufficiently high degree of anisotropy Jï^| , and the number of ions 
in the trap is not too large, the ions will lie in a string along the x-axis. The vibrational motion of the different ions 
will be strongly coupled due to the Coulomb repulsion, and so their motion, when small, is most easily described 
in terms of a set of normal modes. The mode angular freqüències are given by ^/Jt^uj x , where \i p is the eigenvalue 
associated with the mode Q , and ui x is the angular frequency characterizing the strength of the trapping potential in 
the x direction. For simplicity, we will only consider oscillation along the x direction, oscillations transverse to this 
direction are important only for low anisotropies if necessary the theory presented here can be extended to this case 
without too much dimculty). The vibrational quantum state of the ions is specified by a set of excitation numbers 
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for the various phonon modes of the oscillation. If there is a total of N ions in the trap, then there will be N such 
modes, which can be enumerated in order of increasing eigenfrequency (see ref. Q for more details). Each of these 
modes can be considered as an uncoupled quantum harmònic oscillator, and so each has a ladder of number states. 
Thus the state will be denoted |{^}), where {n}ni, nz, ■ ■ ■ is the set of mode excitation numbers, i.e. there are n\ 
quanta in the first mode, n% quanta in the second mode, and so on. 

We will use a single index to denote the combined wavefunction for the internal state and for the vibrational 
state. Thus \i) will stand for I7, J,mj) <8> |ni,ri2, ■ ■ ■ tin). We shall use {\i), \j), \ f) etc.} to denote states in the 
lower manifold of atòmic states and {\k), \ l), \ k'}, \ l') etc.} to denote states in the upper manifold of atòmic states. 
The energies of these states will be denoted Ticüí, hcük etc. Thus the state function for the ion will be given in the 
Schròdinger picture by the following expression 

IV>(*)> = '52a j (t)exp{-L·j j t)\j) +J2 b k(t)exp(-iu k t)\k), (II.l) 
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where a,j is the amplitude in a lower manifold state \j), and bk is the amplitude in an upper manifold state |fc), 
the summations being over all lower manifold states and all upper manifold states, respectively. We will make the 
'essential states' approximation, i.e. assume that these states are the only ones which will have a non-zero amplitude; 
thus the following equations apply: 

Eb')0'i+Ei fc )( fc i| ( IL2 ) 

H Q = ?i|ç^b·)<. ? |+Ç^.|fc)(fc||, (II.3) 

where / is the identity operator and Hq is the unperturbed Hamiltonian operator. 

The interaction Hamiltonian describing the interaction of the atòmic system with the làser field is given by 

Èj = -d-E(f,í), (II.4) 

where d is the dipole moment operator for the ion (which acts only on the internal atòmic states), E is the elèctric 
field strength, which is a function of time t and of the position of the ion, which is expressed in terms of its operator 
f , which acts on the vibrational states. The elèctric field will be assumed to consist of a sum of monochromatic 
components, 

E(r,í) = iïe|Ee (a) S( a )(r)exp(iw a í)| , (II.5) 



where e( a ) is the unit polarisation vector for the a — th làser beam. Substituting eq. (11.5) into eq. (II. 4), and 



using (II. 2) and the fact that, by assumption (fc|d|/s') = (j|d|j') = 0, and making the rotating wave approximation 
(i.e. neglecting terms oscillating at freqüències LUk — luj + uj a ) one obtains the following formula for the interaction 



Hamiltonian: 



ff/ = -EE S n ií b'X fe l eM^J) + h.a., (11.6) 
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where 

- (j\d-e^E^{r)\k). (II. 7) 
Thus one obtains the following differential equations for the amplitudes üj (t) and bk (t) : 

<%•(*) - ^EE^S^pÍ-^W*)' ( IL8 ) 

k a 

= ^EE^tex P (iAgí) a/ (í), (II.9) 



j> a 
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where the asterisk denotes complex conjugation and 

A$ = u k - Uj - u a , (11.10) 

(see fig.2.1). 

If we adiabatically eliminate the populations in the upper manifold (see Appendix A), one obtains the following 
formula for the amplitudes in the lower manifold: 

à l (t) = iJ2^(t)a J (t), (11.11) 

3 

where the coupling matrix /3,j (í) is given by the following formula: 

M) = E E l \(a-) eX P P K ~ *] ' ( IL12 ) 

k a,a' 4A jfe 



B. Raman transitions between two levels with two làser flelds 



We will now examine the special case in which there are only two làser beams present, the "pump" fi eld an d the 
"Stokes" field, with freqüències uj p and uj s respectively, such that uj ps =u) p — oj s > 0. Thus using eq. (11.12), and 
neglecting "rapidly" oscillating terms, (i.e those proportional to exp [±i (u) ps ) t]), one obtains the following expressions 
for the elements of the matrix j3ij (t) : 
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exp [ií(w ps 



(Wj > Ui), 



exp [-ií(w ps - |wíj|)] (uJi > Uj). 



(11.13) 



The constant diagonal elements represent an "A.C. Stark" shift, i.e. a change in the resonance frequency of a transition 
due to an applied optical field. The off-diagonal elements represent induced transitions between pairs of levels of the 
lower manifold. 

We will now make a further specialisation by considering only two lower manifold states. We will assume that the 
quantity lü ps — is large for all but a single pair of states which we will denote \I) and \ J), where (wj > u>i). 
Thus the off-diagonal elements of the matrix will be very rapidly oscillating except for these two states, and there 
will be negligible probability of any transitions being induced. We will denote the average detuning of the fc-th upper 
manifold state by A& = uju — (u>i + wj + lü s + u> p )/2, and the Raman detuning by S = oj ps — ujji (see fig.2.2). If one 
makes the assumption that A& 5 and Aj, 3> u) ps , then the elements of the 2x2 matrix fiu can be written as: 
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The amplitudes for the two states \I) and | J) therefore obey the following equation 



d_ 

dt 



ai 
aj 



A Bé ts 
B*e~ itS D 



aj 



These are equivalent to the well-known Bloch equations for a two level system ||, [ïï 
resonance condition 

8 = A-D. 



(11.17) 

If one makes the following 
(11.18) 



which can be a chieved in practice by varying the relative frequency of the two làser beams, then the solution of 
equation ( 11.17 ) is equivalent to the Standard on-resonance solution of the Bloch equations. Using this solution wc 
can obtain the following solutions for the amplitudes aj(í) and aj(t): 



ai(t) 
aj(t) 



exp 



MA + DV/21 í eÍXcos0 / 2 ie i0 sin0/2 \ / oj(0) 



where 



2|n e// |í 
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and the effective Rabi frequency for Raman transitions between levels |J) and | J) is given by the formula 



fi 



eff 



E 



4A fc 



(11.19) 

(11.20) 
(11.21) 

(11.22) 
(11.23) 



Thus by controlling the duration of each làser pulse precisely, one can control the populations in the two levels \I) 
and | J). Furthcrmore, by controlling the relative phases of the pump and Stokes làsers one can control the relative 
phases of the two states. The phase change x, caused by the A.C. Stark shift, will be small, but must be taken into 
account if large scale quantum computations are ever to be attempted in ion traps. 

The value of fi e / / can be calculated in terms of atòmic constants as shown in Appendix B. In quantum computation 
the phonon modes are assumed to be unpopulated. Two types of transition are of particular use, namely the "V" 
type or carrier transition, for which phonon state is left unchanged and the "U" type or red sideband transition, in 
which one phonon is created M. To be more specific, for a "V" pulse, 



\I) = |0) <g> \vac) 
\J) = |1) <g> \vac), 



(11.24) 



where |0) and |1) are the two infernal levels of the ion which form the qubit, and \vac) is the phonon vacuum state 
(i.e. every mode has zero occupation number); for a "U" pulse, 



V) 



|0)®|{1,0, 
|1) ® \vac), 



,0}) 



(11.25) 



where |{1, 0, . . . , 0}) is the state in which the first phonon mode has occupati on nu mber 1, and all the others have 
occupation number zero. Pulses of these two types for which 6>,defmed in eq. (II. 2C), takes the vales 7r/2,7r and 2tt 
can be combined to constitute a quantum lògic gate jjj. 



III. COOLING, QUANTUM LOGIC AND READOUT BY RAMAN TRANSITIONS 



We will now discuss how the various làser interactions described theoretically in the previous section can used to 
perform the operations required for quantum computing. These operations can be divided into four parts, namely 
Doppler cooling of the trapped ions to produce a crystal (this does not, in fact, require Raman transitions); sub- 
Doppler cooling to reduce the crystal to the quantum ground state of its motion, thereby preparing the quantum 
register in the required initial state; quantum lògic operations; and read-out of final quantum state. 
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A. Doppler cooling 



Use of Standard Doppler cooling [[Ï7| is necessary in order to crystallise the trapped ions into the the required string 
configuration with a temperature small enough to allow sub-Doppler cooling. Doppler cooling is done as follows: 
photons from a làser at 397 nm, detuned to a slightly longer wavelength (i.e. to the red) of the 4 2 S;iy 2 to 4 2 P 1 / 2 
transition are absorbed by the ion, whose momentum is reduced by the conseqüent recoil; this process is followed by 
spontaneous emission, which causes the ion to recoil in a random direction. Thus the wholc process has the net effect 
of, on average, reducing the momentum of the ion along the direction of the cooling beam. 

Because the upper state of the 397 nm resonance, the 4 2 P 1 / 2 state, can also decay to the long lived 3 2 D 3 / 2 state, 
another làser at 866 nm must also be employed to pump out any population which might get trapped there. (see Fig. 
3.1). 

The lowest temperature attained via Doppler cooling, which is limited by the random recoil of the ions during 
spontaneous emission, is given by Tu op = Ti/2KbTps (where Tps is the radiative lifetime of the cooling transition 



and Kb is Boltzmann's constant) |18 . For Ca + , Toop = 4.96 x 10~ 4 K (which, for a trapping frequeney of lj x — 
(27r)x500 kHz, corresponds to an average center-of-mass phonon number of n = KsToop/^x — 20.7). At this 
temperature the trapped ions are crystallised (crystallisation can and does oceur above the Doppler limit). 



B. Sub-Doppler Sideband Cooling via Raman Transitions 



Once a string of ions has been cooled to the Doppler limit, a different cooling technique must be employed to 
remove the remaining phonon excitations of the string. One method of doing this is to perform a U-type red sideband 
transition between the 4 2 Si/ 2 state and the 3 2 D 3 / 2 state using a 732 nm làser pulse detuned by the trapping frequeney 
lü x ; population in the 4 2 S;iy 2 ground state would be excited to the 3 2 D 3 / 2 excited state, and, at the same time, one 
phonon excitation would be annihilated. A V-type carrier transition at 866 nm could then excite the ion to the 
4 2 Pi/ 2 state, from whence it would rapidly decay to the 4 2 Si/ 2 ground state. Provided the trapping potential is 
sufhciently strong, there will be only a small probability of the decay process re-exciting the phonon by recoil (this is 
the Lamb-Dicke effect, discussed below). Thus the ion would be returned to its initial state, allowing the process to 
be repeated until all of the phonons are pumped out. This is a simplified description of single làser sideband cooling, 
which does not take into account the Zeeman structure of the various levels. More details can be found in ref . Q| . 

This method, however, requires high frequeney stability of the làser in order to resolve the motional sidebands. Using 
Raman transitions for the sideband cooling, the làser stability requirement is dramatically reduced |Ï2] |. Instead of 
generating a 732 nm làser with a bandwidth of a few kHz, all that is required in the Raman case is that the relative 
frequeney difference between the Pump and Stokes beams be stable. Sideband cooling using this method consists of 
two steps. First, a "U" type Raman transition is carefully tuned to resolve motional sidebands, allowing transitions 
between the two sublevels of the ground state to be performed while at the same time reducing the phonon number 
by resolving the appropriate sideband. Fig. 3.2a shows this transition as proceeding from the mj = +1/2 to the 
mj = —1/2 sublevels of the ground state. If there was any population in the mj — —1/2 at the start, then this 
process would have the opposite effect, i.e. it would add a quantum of motion rather than subtracting one. For this 
reason we must also carry out the second step, fig. 3.2b, which involves optically pumping the ion back to the original 
sublevel of the ground state. This is achieved by using on-resonance circularly polarised light: any population in the 
4 2 Si/ 2 , m,j = —1/2 sublevel will get excited to the 4 2 P]y 2 , mj = +1/2 sublevel, from where it can decay either to 
the 4 2 S!/ 2 mj = —1/2 or mj — +1/2 sublevels. In the former case, it will get re-excited back to 4 2 P 1 / 2 , mj = +1/2, 
in the later case it will remain, as the circularly polarised light does not have any effect on the 4 2 S 1( / 2 , mj = +1/2 
population. After 100 nsec or so (i.e. ~ 10 cycles), over 98% of the population will end up in the 4 2 Si/ 2 , mj = +1/2 
sublevel, and so step 1 can be repeated. In order to maintain a high cooling rate the 3 2 D 3 / 2 state must again be 
pumped out as it has a chance of becoming populated during the optical pumping process; this will require modulation 
of the 866 nm làser so that all three of the Zeeman levels of the 3 2 D 3 / 2 can be re-pumped. 

The Lamb-Dicke effect will stop the ions from being re-heated during the pump-out process (which involves spon- 
taneous emission and therefore the possibility of heating due to recoil). Specifically, if the kinetic energy imparted 
to the ion by the spontaneous emission recoil is less than the energy required to excite a phonon, then the recoil is 
absorbed by the ion trap apparatus as a whole rather than by the individual ion. Hence, for low phonon numbers, we 
require that 

^ « hu» (HI.1) 
where p = fïk is the momentum due to the recoil (k being the wavevector of the emitted photon), M is the mass of 
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the ion stri ng (6 .64 x 10 26 kg per ion for Calcium), and ui x is the angular oscillation frequency of the ion trap. The 
incquality (III. 1) can be re-written as follows: 



where r\ is the Lamb-Dicke parameter, a dimensionless parameter characterizing the coupling between the làser and 
the motional degrees of freedom of the ion. The value of r\ depends on the transition involved; its effective value for 
Raman transitions is discussed in Appendix B. This simplificd cxplanation using classical mechanics is, of course, 
not vàlid for the quantum case; a thorough analysis shows that the probability of a recoil phonon being created 
during the spontaneous emission is approximately equal to the square of Lamb-Dicke parameter, nrj 2 , where n is the 
average phonon number. Hence for small vàlues of rj and low phonon numbers, the probabilit y of rehcating will be 



small. Since |k| = 2n/X, where A is the wavelength of the emitted photon (i.e. 397 nm), eq. ( III. 1 implies that the 
trapping frequency lo x must satisfy the inequality lu x S> n(27r)31.6[kHz]/N (where N is the number of ions), in order 
for sideband cooling to be efhcient. 



C. Quantum Logic 

Quantum lògic operations can be performed using Raman transitions. As mentioned above, two types of operation 
are required: V-type carrier transitions (fig.3.3(a)) and U-type red sideband transitions (fig. 3.3(b)). Selection 
between U-type and V-type pulses can be made by varying the relative detuning between the Stokes beam and the 
Pump beam. The manner in which these pulses can be combined to realise quantum lògic gates is described in detail 
by Cirac and Zoller in their seminal paper Q, and there is no need to repeat it here. 

As mentioned above, the |1) state has a very low probability of spontaneous emission to the |0) state. Furthermore, 
because there is only a very small probability of exciting the 4 2 Pi/2 level, there is only a very small chance of decays 
to the 3 2 D 3 / 2 level, and so a pump-out làser at 866 nm is not required during quantum lògic operations. 

One problem in performing quantum lògic operations in this manner is the auxiliary state \aux) required to perform 
the Cirac-Zoller CNOT gate. As mentioned in the introduction, 4 2 S!/ 2 the ground state of 40 Ca + has only two 
sublevels, which are employed as the logical |0) and |1) states; a third level is required for \aux). Four possibilitics 
exist. The first is to use one of the sublevels of the other states [^0j: if we use a P state sublevel, there is a danger of 
spontaneous emission while the auxiliary rotation is being performed; if we use a D state sublevel, the quantum lògic 
scheme will require either phase-locked làsers 397 nm and 866 nm to perform Raman transitions between the S to D 
levels, via the P level, or the construction of a narrow bandwidth 732 nm làser. 

The second possibility is to use one of the other sidebands (i.e. a sideband not associated with the CM mode but 
one associated with the breathing mode, for example) [ jÏ3| . However, if the sidebands are close together in frequency 
space that this method will not be particularly easy 

The third way is to use an isotope of calcium which has hyperfine structure. This requires an isotope with non-zero 
nuclear spin, the only stable one being 43 Ca (0.135% natural abundance) with nuclear spin 7/2. Alternatively radio- 
isotopes can be used (see Table 1.), of which 41 Ca or 45 Ca appear to be the most promising candidates. Although 
these isotopes are very scarce, and therefore expensive to purify, it should be noted that ion traps lend themselves 
rather well to use as mass hlters, and so can be used for in situ isotope separation. DeVoe and Kurtsiefer at IBM 
Almadén are performing such isotope purification in their Ba + experiment, although the natural abundance of the 
desired isotope of barium is considerably larger than that of calcium [^T| . 

The fourth possibility is to change the quantum computational paradigm so that an auxiliary level is no longer 
needed. Wineland's group at NIST Boulder has described such an alternative to Cirac and Zoller's quantum lògic 
scheme which has the added advantage of requiring fewer pulses [ p2[ . However this scheme does require precise 
control of làser beam directions and trap freqüències, so that exact and rather large vàlues of the effective Lamb-Dicke 
parameter are attained. 



D. Readout 



An efHcient readout scheme is the one aspect of quantum computation in which ion traps are far superior to any of 
the other proposed technologies. The bàsic idea is to use a "eyeling transition" in which population from one of the 
two quantum levels of each qubit, |1) say, is repeatedly excited by a làser to some short-lived state. This short-lived 
excited state then decays back to the original state |1), resulting in observable fluorescence, without disturbing any 
of the population in the other quantum state |0). This technique was used successfully in the breakthrough quantum 
gate experiment performed at Boulder and reported in 1995 0. 
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In the scheme under discussion here, the readout can be achieved by a two-step process as follows. First any 
population in the 4 2 S!/ 2 , mj = +1/2 sublevel (i.e. the |1) state of each qubit) is transferred to the 3 2 D 5 / 2 level 
by a pulse of a 729 nm làser (Fig. 3.4(a)). This làser does not need to resolve individual motion sidebands of this 
transition, and so it does not have a particularly stringent bandwidth requiremcnt. We shall refer to this as the 
"shelving transition" . 

Secondly, the 4 2 Si/2, mj = — 1/2 to 4 2 P;iy 2 , mj = —1/2 transition is excited using the 397 nm làser (Fig. 3.4(b)). If 
there is any population in the 4 2 S 1 / 2 , mj = —1/2 state (i.e. the |0) statc of the qubit), there will be fiuorescence which 
we can observe. Note that there will be decay from 4 2 P!/ 2 , mj = —1/2 to the 3 2 D 3 / 2 , mj = —3/2, — 1/2& + 1/2 
states, which will trap the populations. Therefore, just as was the case during Doppler cooling, these states must be 
pumped out using 866 nm light. 

If one were using 43 Ca + ions, one might consider an alternative scheme based directly on that successfully used 
by the Boulder group ^]. The 4 2 Si/ 2j F — 3,mp = + 3 and the 4 2 Si/ 2 , F = 4, mp = +4 hypcrfine sublevels would 
form the |0) and the |1) states, respectively. The cycling transition would thcn consist of right-circularly polariscd 
light tuned to excite population in |1) to the 4 2 P 3 / 2 , F = 5, mp = +5 statc. If population excited to this level were 
to decay directly to the ground state, then, because of the electric-dipole selection rules, it could only return to the 
|1) sublevel, resulting in fiuorescence which could be observed. However, the 4 2 P 3 / 2 level can also decay to both 
the 3 2 D 3 / 2 and the 3 2 D 5 / 2 metastable levels. Since re-pumping population trapped in these levels in a manner that 
avoids any false signals is a very complicated problem, this readout scheme can be effectively ruled out. 

Another alternative has been suggested by Steane and his collaborators at Oxford [Q. Instead of performing 
the shelving transition, they propose to optically pump the population of the 4 2 Si/ 2 , mj — +1/2 state via the 
4 2 P 3 / 2 , tíij = +3/2 state to the 3 2 D 5 / 2 level, followed by resonance fiuorescence on the 4 2 S!/ 2 to 4 2 P!/ 2 397 nm 
transition (with the usual 866 nm làser to avoid trapping population in the 3 2 D 3 / 2 level). However the 4 2 P 3 / 2 , mj — 
+3/2 state can decay to both the the 3 2 D 5 / 2 and the 3 2 D 3 / 2 levels, so that the optical pumping scheme has a small 
(10%) chance of causing a false signal. Although this is a drawback, this scheme does have the advantage that a làser 
in the 729-732 nm range is not required. 

Once the readout has been accomplished, the population shelved in the 3 2 D 5 / 2 level has to be removed. This can 
be achieved either by repeating the shelving pulse with the 729 nm làser, or by using two làsers at 866 nm and 854 nm 
to pump population from the metastable state D-states to the ground state via the 4 2 P 3/2 / 4 2 P 1/2 levels. The second 
technique is quicker and more reliable, but it does require a third làser wavelength at 854 nm. 



IV. EXPERIMENTAL REQUIREMENTS 



In this section we will discuss the specifications of the apparatus needed to carry out the experiments described 
above. A full description of the ion trap, the vacuum system, the associated làsers and optical systems can be found 
in ref. ||. Here we will confine ourselves to a description of the additional apparatus needed specifically for the 
stimulated Raman transitions described above. 



A. Static magnètic field 



A strong static magnètic field is necessary in order to split the Zeeman sub-levels in frequency space, so that thcy 
may be reliably resolved with the available làsers. The Zeeman structure of the 397 nm 4 2 Si/ 2 to 4 2 P;iy 2 transition 
and of the 866 nm 4 2 D 3 / 2 to 4 2 P]y 2 transitions are shown in figs. 4.1(a) and 4.1(b), respectively. The minimum 
splittings between component lines are (27r)0.93MHz/Gauss and (27r)0.19MHz/Gauss respectively, while the radiative 
line widths are 129MHz and 10.6MHz, respectively. Thus a magnètic field strength of 200 Gauss will be sufficient for 
the Zeeman splitting of both lines to be well resolved. 

The stability requirements for the static magnètic field are important. Variations in the magnètic field during 
quantum operations will mean that there is an error in the effective Rabi frequency and so there will be an over- or 
under-rotation during the lògic operations. The error in angle of rotation is 




(IV.3) 



Thus if we require a angular accuracy of say 59/9 ~ 10~ 6 (corresponding to an error ratc below the accuracy threshold 
for fault-tolerant quantum computing [p3|) with each operation taking ~10 ^tsec, and we have a magnètic field of 
200 Gauss, then the rate of change of the field must be less than 20 Gauss per second. The magnètic field may be 
stabilised by placing a Hall probe in the field and using it as a sensor for a feedback loop to the Helmholtz coils. 
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B. Làsers 



1. Power 



The power requirements of the difïerent làser tasks needed for quantum computing are approximately given by the 
following formulas (for derivation, see refs. [| and also see Appendix B for Raman transitions): 



( ü 2 hcd 2 T 

riAhcd 2 T 

A 3 



P 



(single làser, V pulse) 
(Raman, V pulse) 



ÜAhcd 2 r ( ] ( Raman; u p U i S e) 



(IV.4) 



Ahcd 1 



(single làser, saturation power). 



In these expressions, Í7 is the effective Rabi frequency for the process, h is Planck's constant, c is the speed of light, d 
is the 1/e 2 diameter of the intensity pattern of the làser in the vicinity of the ion (we have assumed a TEMoo Gaussian 
beam), A and r are the wavelength and radiative lifetime of the transition, rcspectively, A is the Raman detuning, 77 
is the Lamb-Dicke parameter and N is the number of ions. 

For purposes of estimating power, we assumc that we wish to obtain Rabi freqüències equal to (27r)1.0 MHz for 
all the operations. We require tight focusing only for the 397 nm làser during quantum lògic operations; the other 
two làsers (866 nm and 729 nm), and the 397nm làser during Doppler cooling and readout, do not need to address 
individual ions; in fact it is more convenient for them to have a broad focal pattern with intensities equal for all ions 
in the trap. Therefore we assume a, d — 10/im focal spot diameter for the 397 nm làser during lògic operations and 
d = 100/im focal spot diameters for the other làsers. The Lamb-Dicke parameter was assumed to be 77 = 0.1, the 
number of ions to be N = 10 and the Raman detuning to be A = 10 GHz. 

The 397 nm làser will be the workhorse of a Raman-based Calcium ion quantum computer. It will be required to 
perform Raman transitions (of both the U and V types) on single ions for quantum lògic, U-type transitions on all of 
the ions (with a broad focal spot size) for sideband cooling, and to be used to saturate the 4 2 S!/ 2 to 4 2 P 1 / 2 transition 
during the Doppler cooling and the readout phases. We calculated the power required for all of these operations and 
found that the U-type Raman transition during sideband cooling required the most power (which is the figure given 
in Table 2.) The 866 nm làser is needed for pump out, acting on all of the ions (and therefore having a broad focal 
spot); its power is the saturation power for that transition. The 729 nm làser is used for a single V-type ir -pulse 
shelving transition during readout. Because the 4 2 S x / 2 to 3 2 D 5 / 2 transition is dipole forbidden, this làser has the 
highest power requirements. 



2. Optical system 



In order to fulfill momentum conservation and quantum selection rules, the Raman and Stokes beams are required 
to have different polarisations and propagation directions, i.e. each ion must be physically addressed by two làser 
beams, each with slightly different frequency, polarisation, and direction. This complicates the ion addressing òptics, 
but not unduly so: Figure 4.2 shows a simple experimental layout for producing appropriate Pump and Stokes beams. 
The arrangement shown has the disadvantage that the effective wavevector k s — k p has a component transverse to 
the weak axis of the trap, and so care must be taken to avoid the creation of transverse phonons when using it. The 
alternative is to use a optical arrangement in which the pump and Stokes beams are no longer co-propagating, which 
requires duplicate optical addressing system for the two beams. However, only one of the pump or Stokes beams 
needs to be tightly focused to address the ions; the other beam can have a broad focal spot, acting on all of the ions 
simultaneously (and resulting in an A.C. Stark shift of the 4 2 S;iy 2 to 4 2 P x / 2 397 nm transition); transitions would 
only occur when the second, tightly focused làser beam is acting on the ion pjj . 
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3. Làser Stability 

The bandwidths of the 397 nm and 866 nm làsers are determined by the strength of the magnètic field: The Zeeman 
splitting of the transitions, discussed above, must be resolvable in frequency by the làsers. There is no need for làsers 
to be narrower than the natural widths of the lines (which we have used as the bandwidths cited in Table 2. 

For the 729 nm làser needed for the shelving pulse during readout, the stability is determined by the accuracy of 
the 7r-pulse. For a on-resonance làser, all of the population in the lower level is transferred to the upper level by the 
7r-pulse. However, if the làser is detuned by an amount A because of its finite bandwidth, there will be a probability 
(A/fi ) 2 of remaining in the |1) state (and consequently causing a false readout signal). Here Slo is the single làser 
Rabi frequency, which may be calculated using formulas given in ref. For fio = (27r)1.0MHz, a detuning of 
A = (27r)100kHz gives a 0.3% chance of this happening; thus a bandwidth of ~ (27r)100kHz should be sufhcient for 
the 729 nm làser. Such bandwidth are now Standard. A more accurate readout can be achieved by either increasing 
the 729 nm làser power (thereby increasing fio) or by decreasing the bandwidth. 

^. Pump out làser 

As mentioned in section 3, the Zeeman splitting of the 866 nm transition (shown in fig. 4.1) will require some 
form of frequency modulation of the pump-out làser at that wavelength. For a 200 Gauss magnètic field, the amount 
of modulation is roughly 280 MHz. One way to achieve this is with two electro-optic modulators (EOMs) in series, 
one generating frequency sidebands at 18.6 MHz (i.e. Aw/15 ), the other generating sidebands at 224 MHz (i.e. 
12Aw/15). 

5. Radio -frequency oscillator 

Frequency differences between the pump and Stokes beams of several GHz can be achieved with commercial radio- 
frequency (RF) sources and electro-optic modulators. The RF signal needs be stable, the stability required being 
directly dependent on the required accuracy of the lògic operations. If we set this as 1 part in 10 6 we need an RF 
bandwidth of about 1 kHz. Thus the phase stability of commercial RF sources is more than adequate for elementary 
quantum lògic experiments. 

Table 2 shows the làser required characteristics calculated on the basis of the considerations discussed above. Note 
that the calculated làser intensity is for the focal region and so does not take into account losses due to various effccts 
in the optical system. 

V. CONCLUSIONS 

This paper has presented in detail a scheme for performing quantum lògic operations with trapped calcium ions using 
stimulated Raman transitions. Speaking from our own experience in establishing an ion trap quantum computing 
experiment ||, the most challenging single experimental problem is the generation of a large number of powerful cw 
làser beams at a variety of wavelengths and with very stringent stability requircmcnts. This is particularly the case 
when one attempts a "single làser" transition scheme, exploiting the Rabi oscillations between two levels to perform 
the necessary operations. It was this difficulty which was the principal motivation for the approach we have considered 
here: the adoption of Raman type làser transitions means that absolute phase stability of làsers is no longer a critical 
issue. However some of the resultant experimental arrangements are more complex. In particular, there are stringent 
requirements for the stability of the static magnètic field. As experimental quantum computing is without a doubt a 
very challenging problem, one must have as many technological options as possible if one hopes to meet with success. 



Acknowledgments. 

The authors wish to thank Rainer Blatt, Ralph DeVoe, Richard Hughes, Christian Kurtsiefer, Christof Naegerl, 
Andrew Steane, Dale Tupa and David Wincland for useful conversations. This work was funded by the National 
Security Agency. 



11 



APPENDIX A: ADIABATIC ELIMINATION OF THE UPPER LEVELS 



The eliminat ion of the upper state populations can be carried out as follows. The formal solution of eq. ( JII. S| ) is: 

h(t) = M0) + ^£E íl$* f a y (f ) cxp (iA$f ) dt'. (A.l) 

If one now makes the Markov approximation (i.e. assume that ay (í) is slowly varying so that it can be factored out 
of the integral) one obtains the following approximate expression for the amplitude in the upper manifold: 



n 



2A 



(o)* 
j'k 

— a 3' 



exp 



j'k 



(A.2) 



We can assume that the initial amplitude in any of the upper manifold of states is zero, i.e. 6fc(0) = 0. Since all 
of the states of the upper manifold have a dipole allowed transition to one or more states of the lower manifold, it 
is reasonable to assume that any res idua l population in the upper manifold will rap idly decay due to spontaneous 
emission. On substituting from eq. (A.2) into the equation of motion for a,j(t), eq. flH-S| ), and neglecting a rapidly 

oscillating term proportional to exp I 

manifold. 



-iA^í), one obtains the formula eq. ( 11.11 ) for the amplitudes in the lower 



APPENDIX B: CÀLCUL ATION OF THE EFFECTIVE RABÍ FREQUENCY IN TERMS OF ATÒMIC 

CONSTANTS 



From eqs. (H.7) and (11.23), we obtain the following formula for the effective Rabi frequency: 



^ _ y (i\d ■ e(rtg(p)(r)|fc)(fc|d ■ e(">»£;M*(f)|j) _ (R1) 

In order to calculate this quantity in terms of atòmic constants etc, it is necessary to return to the representation of 
the states in terms of tensor products of internal states with external (vibrational) states, vis: 

|») - |0>®|{m}) 
li) - 

\k) -> (B.2) 

where |0) and |1) are two internal states of the lower manifold (i.e. the two states that store the information of the 
qubit of the quantum computer), |A) denotes an internal state in the upper manifold, and |{m}), \{n}} and \{£}) 
denote the state of the external degrees of freedom of the ion (i.e. the vibrational state) specified by a set of occupation 
numbers for the different modes, as discussed above. Note that we have used a simplified notation for the internal 
states, the single symbol (0, 1, or A) denoting the set of quantum numbers 7, J, mj needed to specify the internal 
state completely. 

Using this notation we obtain the following formula for the effective Rabi frequency: 



X i 4 " ( A A - W £ ) 



|d • eW]A)(Ald • e( s )*|l)({m}|£;(p)(f)|{^})({í}|í;( s )*(f)|{r ï }) 

where A\ = lú\ — (lü +lji+lü p +üj s )/2 and = 0J x [^/W{h~ (mi+ni)/2} + sfïíi{li — (m2+n 2 )/2} + . . . ^/JÏn{£n — 
(rriN + n A r )/2}]· We will make the approximation that ^> io^ h \ which is justified in most experimental circum- 
stances, where the detuning A^ is of the order of GHz and the trap frequency oj x is usually is no higher than a few 
MHz. Further, we will assume that both the làser beams are, to a good approximation, plane waves, i.e. 

fiW(f) = E { P) exp (-ik p ■ f ) 

E^ s) (f) = E { s) exp (-ik s • f ) . (B.4) 
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Thus, if we use the completeness of the vibrational states J^e lí^lKÍ^II = then we get the following expression for 

n.„ = £ W* • ° W • ''•"I" E m ({M) , {„}, ,) , (B.5) 



where 



(m p | exp[i^ p (ój + (B.6) 



In eq. ( |B.(^ ), represents the state of the p-th oscillator mode with occupation number n p , à p (àj > ) is the annihilation 

(creation) operator for that mode and, for the s-th ion in the chain, Ç p = ribi^ / nj,^ 4 , bi P ^ being the s-th element of the 
p-th mode normalised eigenvector (ref. ||, table 2) and, as mentioned above, [i v is the eigenvalue. The Lamb-Dicke 
parameter 77, which is a dimensionless number that characterises the coupling between the làser and the vibrational 
degrees of freedom of the ion, is given by the formula: 



^-^W™^' (R7) 

M being the mass of one of the ions. Note that for stimulated Raman transitions, the Lamb-Dicke parameter depends 
on the projection of the difference of the two laser's wavevectors onto the x-axis. Thus by carefully selecting the 
directions of làser propagation, one has i n pr inciple a large degree of control over this coupling. 

The matrix elements appearing in eq. (Bii) can be calculated in closed form (see: ref. |17| , section V.A). The result 



is 



where 



(m| exp^t + à)]|„) = A^^^^j}^ „ } (e 2 )exp(e 2 /2)exp(i$ m ,„), (B. 



7r[n + min{m, n} + sgn(Ç)(m — n)/2]. (B.9) 



In eqs.( |B.8[ ) and flB.9| ), L®{x) is the Laguerre Polynomial p4j ], sgn(x) is the sign of x, and min{m, rí}(max{m, n}) is 
the smaller (larger) of the two integers m, n. 

For quantum computing we are particularly interested in performing two types of operation, namely the "V" type 
or carrier pulse, for which {n} — {m} and the "U" type or red sideband pulse, for which {n} — 1,0, ...,0 and 
{to} = 0, 0, . . . , 0. For these two ideal cases, one can show that 

trs xrx \ /l + 0(ry 2 ) (V pulse) 

/ ({to}, {n}, r,) = { v/VW+ 0{r]2) (u pulse) (B.10) 



The atòmic matrix elements appearing in eq. ( B.5 ) can also be calculated. The result is: 



Hd-eW|A) = \ r Ai \ J \t l) E í JV 1 ^ ]^-e {a) > (B-11) 
x 1 V Acak 3 ^ \ -m v q m x I 9 y 1 

q— — l 

where c is the speed of light, a is the fine structure constant, the term containing six symbols in round brackets 
is the Wigner 3j symbol |25| , v stands for or 1, a stands for p or s, k is the wavenumber of the transition (i.e. 
k = 27r/A, where A is the wavelength) and A the spontaneous decay rate between the upper and lower manifolds 
(i.e the reciprocal of the radiative decay lifetime) and the vectors e q are the usual spherical basis vectors: ei = 

(-l,i,0)/v^,e = (0,0,l),e_ 1 = (l,i,0)/V2. 

Thus, finally we obtain the following formula for the effective Rabi frequency of Raman transitions: 

Ah cak 6 í -r / A\ 
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where (3\ is a dimensionless factor of order unity dependent on the quantum numbers of the upper and lower states 
and on the polarisations of the two làsers: 



0x 



q,q'=-l 



Jo 1 Ja 

-Mq q M\ 



Ji 1 Jx 
-Mi q' M\ 



(B.13) 



In order to estimate power requirements, we assume that the pump and the Stokes làsers are of roughly the same 

2 



intensity, so that E^* 



E, 



(p) 



= 64Pah/e 2 cP , where P is the làser power and d is the 1/e 2 diameter of the 

làser (we have used rcf. |26|], eq. (14.5.27) to derive this formula). Thus if we assume that the Raman detuning is much 
larger than the Zeeman splitting of levels in the upper manifold, so that A\ w A, we obtain 



n 



PA 3 



e// hcrd 2 A 



({m},{n},ri), 



(B.14) 



where T is a factor of order 1. Using this formula, we obtain the formula for làser power used in section 4.2. 



FIGURES 




Figure 1.1. Low-lying lcvels of the singly ionised calcium ion, showing all of thc Zccman-split magnètic sublcvcls, 
the absence of hyperfine structure. The spacing is not to scale. 
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Figure 2.1. Notation used in analysis of Raman transitions in the general case. 
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Figure 2.2. Notation used in analysis of Raman transitions in two-level systems. 
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Figure 3.1. Doppler cooling. 
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Figure 3.2(a). Raman sideband cooling, step 1: removal of one phonon. 
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Figure 3.4(b). Raman readout, step 2: resonance fluorescence at 397 nm. 
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Figure 4.1. The Zeeman splitting of (a) the 397 nm 4 2 Si/ 2 to 4 2 Pi/ 2 transition and (b) of the 866 nm 4 2 D 3 / 2 
to 4 2 P 1 / 2 · The quantity Au; is given by the formula /j-bB/H, where [Ib is the Bohr magnetron; numerically, Aw = 
(27r)1.4[MHz]B[Gauss]. The resonances represent the following transitions (from lcft to right, quantum number of 
lower sublevel stated first): (a) mj = +1/2 to mj = -1/2; +1/2 to +1/2; -1/2 to -1/2 and -1/2 to +1/2; (b) 
mj = +3/2 to mj = +1/2 ; +1/2 to -1/2; +1/2 to +1/2; -1/2 to -1/2; -1/2 to +1/2 and -3/2 to -1/2. 
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Figure 4.2. Simple experimental layout for producing appropriate Pump and Stokes beams. 
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TABLES 

Table 1. Isotopes of Calcium with non-zero nuclear spin and half-lives r greater than 1 second. 



Isotopc 


Decay Mode 


Half Life, r 


Nuclear Spin, / 


41 Ca 


Electron Capture 


103,000 years 


7/2 


43 Ca 


(stable) 




7/2 


45 Ca 


f3~ emitter 


162 days 


7/2 


47 Ca 


(3~ emitter 


4.5 days 


7/2 


49 Ca 


/3~ emitter 


8.72 min 


3/2 


51 Ca 


/3~ emitter 


10 sec 


3/2 



Table 2. Estimated làser wavelength, power and bandwidth requirements for performing thc operations described in 
this paper. Numbers are intended for general guidance only. See text for explanation of assumptions and parameters 
used. 



Làser 


Power 


Bandwidth 


397 nm 
866 nm 
729 nm 


0.5 mW 
0.1 a*W 
0.2 W 


(2tt) 20 MHz 
(2?r) 2.0 MHz 
(2tt) 100 kHz 
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